Supersonic strain front driven by a dense electron-hole plasma 
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We study coherent strain in (001) Ge generated by an ultrafast laser-initiated high density 
electron-hole plasma. The resultant coherent pulse is probed by time-resolved x-ray diffraction 
through changes in the anomalous transmission. The acoustic pulse front is driven by ambipolar 
diffusion of the electron-hole plasma and propagates into the crystal at supersonic speeds. Simu- 
lations of the strain including electron-phonon coupling, modified by carrier diffusion and Auger 
recombination, are in good agreement with the observed dynamics. 

PACS numbers: 61.10.Nz, 63.20.-e, 42. 65. RE 
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Subpicosecond laser-induced electron-hole plasmas in 
semiconductors can produce large amplitude lattice 
strain and rapid loss of translational order. These effects 
have been studied extensively in ultrafast linear and non- 
linear reflectivity experiments |[ ||, |[ [| [?|, |], |t| and, 
more recently, in time-resolved x-ray Bragg scattering 
experiments 0, |T| [t|, X-ray diffraction has 

the advantage that it can provide quantitative structural 
information. Many of the x-ray experiments jioj 11 
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have been analyzed using the thermoelastic model put 
forward by Thomsen et al. JtJ in which the strain is caused 
by differential thermal expansion. Deviations from this 
model are discussed in the work of Thomsen et al. and 
have been seen in x-ray diffraction |ll], [l| [l6| fl7| . 
Cavalieri et aZ.|l^, Q studied coherent strain near the 
thermal melting threshold in Ge and concluded that the 
strain is produced over a region which is thick compared 
to the optical penetration depth due to ambipolar dif- 
fusion. However, their experiment was only sensitive to 
structural changes in the near surface region. 

In this letter we report on measurements of coher- 
ent strain generation in Ge following ultrafast laser- 
excitation using a bulk sensitive structural probe. We 
use time-resolved ultrafast x-ray transmission to measure 
strain propagation deep within the crystal, providing in- 
formation about the generation process. Initially, the 
strain front advances at speeds greater than the sound 
speed. In our experiments, the laser intensity is sufficient 
to impulsively generate a dense electron-hole plasma at 
the crystal surface, the dynamics of which are governed 
by ambipolar diffusion Jl8[ and Auger recombination. 
The plasma couples to the lattice through the deforma- 
tion potential. In order to probe the resulting coherent 
acoustic pulse as it travels deep within the bulk, we uti- 
lize the Laue geometry whereby the x-rays traverse the 
full thickness of the crystal, emerging on the other side 
as two mutually coherent beams |Q9[. We have recently 



shown that a short acoustic pulse can coherently transfer 
energy between these two beams on a time scale inconsis- 
tent with the thermo-elastic model. Following the initial 
transient, the beam intensities oscillate as a function of 
the pump-probe delay In the experiments reported 
here, the strain generation is studied as a function of the 
incident laser fluence. The relative phase of the oscilla- 
tions and the amplitude of the transient provide informa- 
tion about the strain generation process at times shorter 
than the x-ray probe duration. 

In the Laue geometry, two linearly independent wave 
solutions propagate through the crystal. Transverse to 
the propagation, these two solutions are standing waves 
whose wavelengths are twice the spacing of the diffracting 
planes. The solutions are usually labelled a and (3 with 
the convention that a has its nodes, and p its antinodes 
on the diffracting planes. In the case that all atoms lie on 
these planes, a is maximally transmitted and (3 is maxi- 
mally absorbed. Because the two solutions interact with 
different electron densities, they propagate with different 
velocities. 

Outside the crystal, two diffracted beams are pro- 
duced: one in the direction of the input beam (forward- 
diffracted or "0" beam), and the other in the direction 
determined by the vector sum kn = fco + Gh (deflected- 
diffracted or "H" beam). Here k^kn) corresponds 
to the wavevector of the forward-diffracted (deflected- 
diffracted) beam and Gh is the reciprocal lattice vector 
corresponding to the diffracting planes. These beams are 
linear combinations of the two internal solutions, a and 
P. The external intensities are given by: 
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where Io(Ih) is the diffracted intensity of the forward 
(deflected) beam, E a ^ is the complex wave field inside 
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FIG. 1: Intensities of the forward- and deflected-diffracted 
beams (upper) and the interior solutions a, (3 (lower) as a 
function of depth inside a thick crystal. A lattice disturbance 
near the exit couples the two solutions, regenerating (3 . 
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FIG. 2: Experimental setup. 



the crystal, k a ^ is the complex wavevector of the a, (3 
solutions (including absorption), and a,b,c,d are deter- 
mined by the crystal orientation. The two internal modes 
E a ,p oscillate in and out of phase as they propagate 
through the crystal. The wavelength of the interference, 
A = \k a — fc/jl -1 , is known as the Pendellosung length 
which is typically a few to tens of microns and is often 
shorter than the absorption length. 

For a crystal that is thick compared to the (3- 
absorption length, only the a solution survives and there 
are no interference effects. This is the anomalous trans- 
mission of x-rays, known as the Borrmann Effect [jl9|| . A 
distortion of the lattice can cause a redistribution of the 
interior wave solutions |2(J . Figure [l] shows the effect of 
a thin region of distortion regenerating the [3 solution af- 
ter it has decayed away in a thick crystal for the case of 
zero a absorption. When this occurs close enough to the 
crystal exit, the regenerated /3 wave does not decay away 
and interference occurs at the exit face, despite the fact 
that the crystal is thick. 

In our experiments, a short acoustic pulse is generated 
at the surface of a thick crystal. This pulse can be con- 
sidered as a moving lattice disturbance. The diffracted 
intensities will oscillate in time as the pulse travels into 
the crystal bulk with a period that is given by the Pen- 
dellosung length divided by the speed of sound. Devia- 
tions from the impulsive strain generation will be evident 
in the phase and/or amplitude of the x-ray modulation 
as a function of pump-probe delay. 

The experiments were performed at the 7-ID undula- 



tor beamline at the Advanced Photon Source. The x-ray 
energy was set to lOkeV using a cryogenically cooled Si 
111 double crystal monochromator leading to a 1.4x 10 -4 
fractional energy spread. The x-ray beam is masked by 
tantalum slits to ensure that the x-ray spot is smaller 
than the laser spot on the sample and to provide x-ray 
collimation. The sample is a 280/im thick, (001) Ge sin- 
gle crystal. The crystal was oriented such that the x-rays 
diffracted in the asymmetric 202 Laue geometry. In this 
geometry, and at 10 keV, the Pendellosung length is 6.2 
/im and the (3 absorption length is 19/im, normal to the 
surface. Therefore, in the unperturbed crystal, only a 
survives at the exit. The only difference between the two 
diffracted beams is in their direction and a mismatch in 
their amplitudes due to details of the boundary condi- 
tions on the exit face [|l9|j . 

Coherent strain pulses are produced on the x-ray out- 
put face of the crystal by sub-lOOfs, 800 nm laser pulses 
at a 1kHz repetition rate. The excitation is fully re- 
versible between subsequent laser pulses. The laser is 
phase-locked to the x-ray beam to better than the x-ray 
pulse duration. The laser is timed to the x-rays using 
a combination of a digital delay generator and an elec- 
tronic phase shifter in the phase-locked loop. In this 
manner the pump-probe delay may be set across a range 
of ±1 ms with 19 ps precision. A fast silicon avalanche 
photodiode (APD) and a picosecond x-ray streak camera 
pl[ were used as the time-resolved detectors. The APD 
sampled the deflected-diffracted beam intensity and the 
streak camera sampled the forward-diffracted beam (see 
Fig. ||). The x-ray bunch separation was ^152 ns, large 
enough to allow electronic gating and measurement of a 
single x-ray pulse. 

Following laser-excitation, high contrast oscillations 
are observed in the pump-probe data over a large span of 
excitation densities. Figure || shows these oscillations in 
the deflected-diffracted beam. The period of oscillation 
agrees with the Pendellosung length divided by the longi- 
tudinal speed of sound. At an incident fluence of 35 JS 3, 
the behavior near t — shows a large transient that is 
unresolved with the 100 ps x-ray probe beam. After the 
transient, the oscillations show a significant phase-shift 
with respect to oscillations that occur following an ex- 
citation of 2^5. The amplitude and frequency of the 
oscillations are relatively insensitive to the fluence. How- 
ever, as shown in Fig. |], the phase is strongly dependent 
on the fluence and is correlated with the amplitude of 
the initial transient. The relative phase of the oscillation 
was defined with respect to the 2^^ excitation and was 
retrieved from a least squares fit |23) . The amplitude 
of the transient is defined as the diffracted intensity at 
a delay of 200 ps. Most of the energy transfer occurs in 
^40ps, measured with the forward diffracted beam us- 
ing a streak camera (see the inset in fig.||). At relatively 
high fluences (> lO- 2 ^), the intensity of the deflected 
diffracted beam approximately doubles while the forward 
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FIG. 3: Time-resolved anomalous transmission. The time- 
dependent intensity of the deflected-diffracted beam at three 
different incident optical fluences: 35 (solid line), 
(dashed line), 2^4^ (dot-dashed line). Inset: Streak cam- 
era data showing the intensity of the forward-diffracted beam 
with picosecond resolution at an incident optical fluence of 35 




FIG. 4: The relative phase of the Pendellosung oscillations 
(squares) and the normalized deflected-diffracted intensity at 
a time delay of 200ps (circles) as a function of incident optical 
fluence. 

diffracted beam is cut by more than 75%. At relatively 
low fluences (< 2-^), there is no transient. 

Inspection of (|l| and (|^) shows that the maximum en- 
ergy transfer between the forward and deflected beams 
near the exit of a thick crystal occurs if the a and 
solutions are coupled at a depth of A/4. This implies 
that the transient behavior is due to a perturbation to 
the lattice that reaches a depth of more than 1.5 /im into 
the bulk. In the simplified picture that a moving inter- 
face couples the a and /3 solutions, the excitation must 
propagate into the bulk at greater than 37,000 m/s, more 



than seven times the longitudinal speed of sound |24j. 

The strain pulse has a finite spatial extent and is com- 
prised of a spectrum of phonons with different wavevec- 
tors. We expect that the phonon component with wave- 
length equal to the Pendellosung length will resonantly 
couple the two interior wave solutions (25) . To model this 
phenomenon, we solve the equations of dynamical diffrac- 
tion within the crystal, taking into account the laser- 
induced time-dependent strain profiles, using the Takagi- 
Taupin formalism adapted for Laue geometry |p6[ |27| . In 
this method, the differential equations coupling the a and 
P branches are solved numerically. The depth-dependent 
strain profile for a given time is taken into account by 
noting that local strain is equivalent to a change in the 
local Bragg angle. Details of this approach (for Bragg ge- 
ometry) can be found in the original work of Takagi|26| 
and Taupin p7[ . The means by which the method can 
be adapted for Laue geometry are implicit in the work of 
Zachariasen |2|] and Batterman and Cole Jl9| . 

Pure thermoelastic models of strain propagation do 
not predict the observed fluence dependence of the phase 
and amplitude of the Pendellosung oscillations. A proper 
model must include the effects of the coupling of the pho- 
toexcited plasma to the crystal lattice. The strain is 
comprised of both diffusive and elastic components: the 
diffusive strain is determined by the instantaneous tem- 
perature and carrier density profiles, modified by Auger 
recombination; the elastic strain is driven by changes in 
the temperature and carrier density, and propagates into 
the crystal at the speed of sound. In the absence of diffu- 
sion, a bipolar pulse develops in the time given by the op- 
tical penetration depth divided by the speed of sound [Q. 
For LA phonons with wavectors along [100] and a 0.2 /im 
penetration depth, this corresponds to ~40ps. Including 
diffusion, the electron-hole plasma extends ~1 /im in the 
same time, leading to a strain front that has propagated 
into the bulk faster than the speed of sound. 

Figure |5| shows the calculated diffraction intensity as 
a function of laser delay at an absorbed laser fluence of 
3^ (corresponding to a carrier density of ~ 6 • 10 20 
cm -3 ). Good qualitative agreement with the experiment 
is seen (Fig. ^). The sharp initial rise in diffraction in- 
tensity is reproduced, as well as the frequency and phase 
of the time-resolved Pendellosung oscillations. Figure ^ 
shows the calculated phase shift and the diffracted in- 
tensity as a function of absorbed optical fluence. After 
taking into account the surface reflectivity of the sample, 
Ijpod agreement with the experiment is obtained (Fig. 

In conclusion we have demonstrated a bulk sensi- 
tive probe of lattice dynamics using time-resolved x- 
ray anomalous transmission. We have observed that 
electron-phonon coupling modified by carrier diffusion 
is a dominant mechanism for energy transport in laser- 
excited Ge. This work could be extended to study how 
the elastic response of the material can modify the elec- 
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FIG. 5: Simulated deflected-diffracted intensity for an ab- 
sorbed laser fluence of 3^4. 
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FIG. 6: The calculated relative phase of the Pendellosung 
oscillations (squares) and the calculated normalized deflected- 
diffracted intensity at a time delay of 200ps (circles) as a 
function of absorbed optical fluence. 



tronic transport properties of semiconductors. 

We thank Bernhard Adams, Marcus Hcrtlcin, Don 
Walko, and Jared Wahlstrand for technical assistance 
and stimulating discussions. We also thank Jin Wang 
for use of the intensified CCD camera. This work was 
conducted at the MHATT-CAT insertion device beam- 
line at the Advanced Photon Source and was supported 
in part by the U.S. Department of Energy, Grants No. 
DE-FG02-99ER45743 and No. DE-FG02-00ER15031, by 



the AFOSR under contract F49620-00- 1-0328 through 
the MURI program and from the NSF FOCUS physics 
frontier center. One of us (SF) acknowledges the financial 
support of Science Foundation Ireland. Use of the Ad- 
vanced Photon Source was supported by the US Depart- 
ment of Energy Basic Energy Sciences, Office of Energy 
Research under Contract No. W-31-109-Eng-38. 



[1] C. Shank and D. Auston, Phys. Rev. Lett. 34, 479 (1975). 
[2] C. V. Shank, R. Yen, and C. Hirlimann, Phys. Rev. Lett 

51, 900 (1982). 
[3] H. W. K. Tom, G. D. Aumiller, and C. H. Brito-Cruz, 

Phys. Rev. Lett. 60, 1438 (1988). 
[4] K. Sokolowski-Tinten, J. Bialkowski, and D. von der 

Linde, Phys. Rev. B 51, 14186 (1995). 
[5] K. Sokolowski-Tinten et al., Phys. Rev. B 58, R11805 

(1998). 

[6] C. Thomsen et al., Phys. Rev. Lett. 53, 989 (1984). 
[7] C. Thomsen et al., Phys. Rev. B 34, 4129 (1986). 
[8] N. V. Chigarev et al., Phys. Rev. B 61, 15837 (2000). 
[9] O. B. Wright et al., Phys. Rev. B 64, 081202 (2001). 
C. Rose-Petruck et al., Nature 398, 310 (1999). 
A. M. Lindenberg et al., Phys. Rev. Lett. 84, 111 (2000). 



[10 
[11 

[12] C. Siders et al., Science 286, 1340 (1999). 
[13] A. Cavalleri et al., Phys. Rev. Lett. 85, 586 (2000). 
[14] A. Cavalleri et al., Phys. Rev. B 63, 193306 (2001). 
[15] K. Solokolowski-Tinten et al., Phys. Rev. Lett 87, 225701 
(2001). 

[16] D. A. Reis et al., Phys. Rev. Lett. 86, 3072 (2001). 
[17] M. F. DeCamp et al., Nature 413, 825 (2001). 
[18] J. Young and H. van Driel, Phys. Rev. B 26, 2147 (1982). 
[19] B. Batterman and H. Cole, Rev. Mod. Phys. 36, 681 
(1964). 

[20] A. Authier, S. Lagomarsino, and B. Tanner, eds., X-ray 
and Neutron Dynamical Diffraction: Theory and Appli- 
cations (NATO ASI Series, 1996). 

[21] Z. Chang et al., Appl. Phys. Lett. 69, 133 (1996). 

[22] H. Barkhudsen et al., J. Mag. Res. 61, 465 (1985). 

[23] F. W. Wise, M. J. Rosker, G. L. Millhauser, and C. L. 
Tang, J. Quant. Elec. 23, 1116 (1987). 

[24] Modest strains can reduce the Pendellosung length by 
more than one order of magnitude due to the change 
in the diffraction condition or crystal structure factor. 
However, the associated fringe visability due to this effect 
will be less than 5 percent. 

[25] I. R. Entin et al., Sov. Phys. Solid State 20, 754 (1978). 

[26] S. Takagi, Acta Cryst. 15, 1311 (1962). 

[27] J. Burgeat and D. Taupin, Acta Crystallogr. A 24, 99 
(1968). 

[28] W. H. Zachariasen, Theory of X-ray Diffraction in Crys- 
tals (John Wiley and Sons, Inc., 1945). 



